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1. The expected phase diagram of QCD 

It is belived that the phase diagram of QCD is very rich, especially at low temperatures, where 
many different phases separated by phase boundaries are expected. The analysis which lead to these 
expectations are based on symmetries, effective theories, models, and analogies, and therefore are 
far from being well established [|Xjl . 

The major problem in the study of the QCD phase diagram is that standard Monte Carlo 
simulations cannot be used in QCD with baryon chemical potential since the fermion determinant 
is complex and cannot be included in a probability measure. It is the so called sign problem. But 
even at zero chemical potential the situation is not clear: in QCD with two degenerate light quarks 
the nature of the deconfinement transition at finite temperature is not definitely established, and it is 
not even clear whether it is a crossover or a true phase transition, either of first or second order class 
[gj. Until recently, it has been generally accepted that it is a crossover, in which case the change 
from the confinement to the deconfinement regime at high temperature and low baryon chemical 
potential would also be a crossover. Since a first order deconfinement transition is expected at 
sufficiently low temperatures, there must be a critical point located where the first order transition 
ceases to exist and is replaced by a crossover. This is the critical endpoint. The expected phase 
diagram is depicted in Figure 1. A Chiral Random Matrix Model predicts this phase structure [Q]. 
Lattice computations overcoming the sign problem by means of the double reweighting technique 
have found signals of a critical point [Q]. 

In the following we will propose a method to locate the critical endpoint and present some 
preliminary results. We assume that at zero chemical potential the confinement and deconfinement 
regimes are separated by a crossover. 



2. The method 



Let us describe the general setting. A more complete discussion can be found elsewhere [9]. 
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Figure 1: The expected phase diagram of two flavor QCD 
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We work on a lattice with a generalized staggered fermion action 

1 3 

S = 9 L L ¥nVi(n) ( U n j\l/ n+i - U^_ u Y„-i I + m£ % Yn + S t , (2. 1) 

where 5 T contains the temporal (co variant) derivative, and therefore the temporal links: 

St = ^X>»I7o(n) (UnfiVn+O-U^oYn-o) + £ £ #,1?o(») feoVn-W) + ^n-O.oVn-o) • (2-2) 

We consider x and y as independent parameters. QCD is recovered by setting x = cosh ft and 
y = sinhjU, where ji is the baryon chemical potential. The sign problem is present if Rey / 0. 
However, it is easy to see that there is no sign problem if y is purely imaginary, y = iy, with y real. 
Imaginary chemical potential corresponds to x = cos/x and y = isin/x. 

It can be seen that the partition function (and therefore the phase diagram) depends on x and y 
only through the following two variables: 

u = x 2 -y 2 , v = (x + y) 3L > + (x-y) 3L > , (2.3) 

where L t is the lattice size in the temporal direction. This two variables are real for both real and 
imaginary y. Also, they are invariant under the change y — > — y (or y — > —y) and x — > — x separately, 
since L t is even. Altogether, this implies that the phase diagram can be plotted in the (x 2 ,y 2 ) 
plane. Obviously, the semiplane y 2 < corresponds to imaginary y. The physical region is the line 
x 2 — y 2 = 1 , with x 2 > 1 . Imaginary chemical potential corresponds to the line x 2 — y 2 = 1 , with 
< x 2 < 1. Numerical simulations are feasible for x 2 > and y 2 < 0. 

The expected phase diagram projected onto the (x 2 ,y 2 ) plane is displayed in Figure 2. The 
extended phase diagram is three dimensional, with the temperature axis, which is not displayed, 
perpendicular to the x 2 and y 2 axis. We assume that on the physical line (the black line on the 
figure, made up of solid and dashed pieces), for large x (i.e, large chemical potential) there is a first 
order phase transition (solid line), which ends at the critical endpoint, marked with an open symbol 
on the figure. For lower x on the physical line there is a crossover (dashed line) separating the low 
temperature confining regime from the high temperature deconfined phase. The first order line is 
extended over a surface in the three dimensional phase diagram, and the critical point is extended 
along a line of critical points (blue line) which is the boundary of the first order surface (projected 
onto the (x 2 ,y 2 ) plane on the figure). 

The y 2 = axis is interesting because numerical simulations are feasible there. At x = 1 and 
y = 0, which is the zero chemical potential point on the physical line, there is a crossover, by as- 
sumption. On the other hand, at x = and y = the fermion determinant contains no temporal link 
-it contains only the spatial links- and obviously cannot break the Polyakov symmetry explicitely. 
Hence, the Polyakov loop is an order parameter. We expect center symmetry spontaneously broken 
at high temperature and a first order transition separating the low and high temperature phases. 
First order transitions are generically robust under perturbations and therefore there must be a first 
order transition for y 2 = and any x 2 sufficiently small. This first order line must end at a critical 
point located at x 2 < 1 (marked with an open symbol on the figure), since we assume a crossover 
at x 2 = 1. The natural scenario is a critical line crossing the physical line at the physical critical 
endpoint and the y 2 = axis at some point of the segment < x 2 < 1, as the blue line of Figure 2. 



161/3 



Searching for the critical endpoint in QCD with two quark flavors 



Victor Laliena 



The critical line can be determined for y 2 < by means of standard Monte Carlo simulations 
and extended to small positive y 2 by double reweighting. The resulting line might be extrapolated to 
the physical line and then the position of the physical critical endpoint would have been determined. 
Even if the extrapolation were not reliable, finding a critical line directed towards the physical line 
would be very interesting, since it would indicate the existence of the critical endpoint. 

3. A Chiral Random Matrix Model 

A Chiral Random Matrix Model of QCD predics the existence of a critical endpoint in the 
(jU,r) plane for light quarks [^]. We can adapt this model to study the structure of the extended 
phase diagram, in the (x 2 ,y 2 ,T) space. To this end, let us consider the following partition function 

3T = J [olAUA\ [dy L dy L dy R dy R ] exp (-S) , (3.1) 

with 



t=0 f=0 



f=0 



z t=0 z f=0 

where A t is a complex N xN random matrix depending on the temporal index t (N collects the color, 
Dirac, and spatial indices), y/£(?) and y/g(t) are respectively the left and right handed components 
of a fermion field, and a is a flavor index. The parameter m is the fermion mass, which will be set to 
zero in the following, x and y are the parameters introduced in the previous section, and j3 controls 
the temperature (large j3 corresponds to high temperature). Notice that, in order to introduce the 
chemical potential in the lattice way, we keep the dynamics local (and free for simplicity) in the 
temporal direction. 




Crossover 



Figure 2: The expected phase diagram projected onto the (x 2 ,y 2 ) plane. Although the transition lines are 
plotted as straight lines for simplicity, it should not be inferred we expect they have such simple form. 
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Figure 3: The physical phase diagram of the CRMM (left) and its extension to the (x 2 ,y 2 ) plane (right). 

The model can be solved analytically in the large N limit with the standard techniques. Fig- 
ure 3 (left panel) displays the physical phase diagram (x = cosh/i, y = sinh/x) for N t = 4. There 
are two phases: at low temperature chiral symmetry is spontaneously broken and it is restored at 
high temperatures. The transition is second order at small pL and first order at large jit. The sepa- 
ration between both lines is a tricritical point. With a small mass, the second order line becomes a 
crossover, the first order transition remains and the tricritical point becomes the critical endpoint. 
The phase diagram mimics perfectly what is expected in QCD. Its extension to the (x 2 ,y 2 ) plane, 
displayed on the right panel of Figure 3, complies with our expectations: the black line is the 
physical line, the solid and dashed red lines are respectively first and second order transitions, and 
the blue line is the tricritical line, which in this case happens to be straight. It crosses the y 2 = 
axis at x 2 = y 2 = 0. This is a pathology of the model which can be easily understood. Moreover, 
the (x 2 = 0, y 2 = 0) point is singular and the tricritical line could not be continued from it by, say, 
double reweighting. Notice however that this model is based only on chiral symmetry and does not 
incorporate the dynamics of the Polyakov loop, which is the basis of the argument of the previous 
section. Hence, it is interesting that its extended phase diagram showed the features we predicted. 

4. Preliminary results for QCD 

Figure 4 displays histograms of the modulus of the Polyakov loop from simulations of QCD 
performed with the R-algorithm. Since it has been pointed out that the systematics of this algorithm 
lead to wrong results for the time steps which use to be used the results should be taken 
with some caution (the extrapolation to zero time step should be addressed before extracting any 
conclusion). Anyway, they are very preliminary. We take two flavors of staggered quarks, with 
mass ma = 0.05, and explore the y = axis on lattices of temporal extent L t = 4 and spatial extent 
L s = 6, 8, 10. For each x the value of the gauge coupling, j3, has been tuned to its "critical" value. 

For x = we see a kind of (not very clear) double pick structure on the L s = 6 lattice, wich 
becomes clearer as L s increases. The data point out to a first order transition (remenber the action 
has exact Polyakov symmetry in this case). The situation for x = 0.5 is similar. For x = 0.7, how- 
ever, it is the opposite. We see a clear two pick structure on the L s = 6 lattice which is disappearing 
as L s increases, signaling a crossover. There must be a critical point between x = 0.5 and x = 0.7. 
Locating it accurately is a hard task which is in progress. 
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Figure 4: Histograms of the modulus of the Polyakov loop in two flavor QCD. 
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